On the motion of three-dimensional compressible isentropic flows 
with large external potential forces and vacuum* 

Jing Li 

f~>) ' Institute of Applied Mathematics, AMSS & Hua Loo-Keng Key Laboratory of Mathematics, 

CN ■ Chinese Academy of Sciences, Beijing 100190, P. R. China 

^ . E-mail: ajingli@gmail.com 

o ■ 

Jianwen Zhang 

School of Mathematical Sciences, Xiamen University, Xiamen 361005, P. R. China 

E-mail: jwzhang@xmu.edu.cn 

Oh ■ Junning Zhao 

■ School of Mathematical Sciences, Xiamen University, Xiamen 361005, P. R. China 

. E-mail: jnzhao@xmu.edu.cn 
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Abstract. We study the global existence and uniqueness of classical solutions to the three- 
dimensional compressible isentropic Navier-Stokes equations with vacuum and external potential 
J> ' forces which could be arbitrarily large provided the initial data is of small energy and the unique 

steady state is strictly away from vacuum. In particular, the solution may have large oscillations 
and contain vacuum states. For the case of discontinuous initial data, we also prove the global 
existence of weak solutions. The large-time behavior of the solution is obtained simultaneously. 
It is worthwhile mentioning that the compatibility condition on the initial data and the regularity 
condition of the external potential forces in the present paper are much weaker than those 
assumed in the existing literature. 
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1 Introduction 

The motion of three-dimensional viscous compressible isentropic flows occupying a domain 
Q C M 3 is governed by the compressible Navier-Stokes equations: 

pt + div(pu) = 0, (1.1) 
(pu) t + div(pn (g> u) + VP(p) = fiAu + (p + A)Vdivu + pf, (1.2) 

where the unknown functions p > 0, u = (u 1 , u 2 , u 3 ) and P(p) = Ap 1 {A > 0,7 > 1) are the 
fluid density, velocity and pressure, respectively, and / = f(x) is the known external force. The 
viscosity coefficients \i and A satisfy the physical restrictions \i > and 3A + 2p > 0. 
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Research Funds for the Central Universities (Grant No. 2010121006), and the Natural Science Foundation of 
Fujian Province of China (Grant No. 2010J05011). 



1 



Let Q = M 3 and poo > be a fixed positive constant. For the external force in the form: 

;=v/ with / = /(*) (1.3) 

we look for the solutions, (p(x,t),u(x,t)), to the Cauchy problem of (jl.ip . (jl.2p with the far 
field behavior 

(p, u)(x, t) — > (pooi 0) as |s| — > oo, t > 0, (1-4) 

and the initial data 

(p,«)(x,0) = (p ,«o)(z), xGM 3 . (1.5) 

The equations (|1.1|) . (|1.2p describe the conservation laws of mass and momentum, respec- 
tively. There has been a lot of literature on the existence and the large-time behavior of solutions 
to the compressible Navier-Stokes equations. The one-dimensional problem has been extensively 
studied, see [HUH] and the references therein. For the multi-dimensional case, the local exis- 
tence and uniqueness of classical solutions was proved in [29^31] in the absence of vacuum, and 
recently, the local strong solutions were studied in [2-4, 30] for the case that the initial density 
need not be positive and may vanish in open sets. The global classical solutions were first ob- 
tained by Matsumura-Nishida [25] when the initial data are close to a non-vacuum equilibrium 
in Sobolev space see also [26J for the exterior problem. Later, Hoff [10|J11| considered the 
weak solutions without vacuum and external forces for the discontinuous initial data. The global 
theory for the multi-dimensional compressible Navier-Stokes equations with "large data" is more 
delicate. Vaigant and Kazhikov [33J studied the global existence of classical solutions in two- 
dimensional periodic domain when the viscosity coefficients depend on density in a very specific 
way and the initial density is strictly positive. One of the most important breakthrough about 
the global theory of "large data" is the work of Lions [24] (see also Feireisl et al. [TIE] ) , who first 
proved the global existence of weak solutions (the so-called "finite energy weak solutions") to 
the initial/inital-boundary value problem of (jl.ip . (|1.2p with generally large initial data when 
the initial energy is finite and the adiabatic exponent 7 is suitably large (i.e. 7 > 3/2). Recently, 
under the additional assumptions that the viscosity coefficients p and A satisfy p > max{4A, — A} 
and that the far field density is away from vacuum, Hoff et al. (cf. [12H14] ) obtained a new type 
of global weak solutions with small energy, which have extra regularity information compared 
with those large weak ones constructed by Lions [24] and Feireisl [HE]. More recently, Huang- 
Li-Xin [T7] established the global existence and uniqueness of classical solutions to the Cauchy 
problem for the isentropic compressible Navier-Stokes equations in three-dimensional space with 
smooth initial data which are of small energy but possibly large oscillations; in particular, the 
initial density is allowed to vanish, even has compact support. 

It has been mentioned in many papers (see, e.g. [91 1221127] ) that the large external forces will 
significantly affect the dynamic motion of compressible flows and cause some serious difficulties 
in the mathematical study. In the following, we briefly recall some recent progress on the multi- 
dimensional compressible Navier-Stokes equations subject to external forces. Indeed, when both 
the initial perturbations and the external forces are sufficiently small, there have been many 
studies on the global existence and the large-time behavior of the compressible Navier-Stokes 
equations, see, for example, [5,6,26,32,34], and among others. For large external forces, Feireisl 
and Petzeltova [9], Novotny and Straskraba [28] proved for different boundary conditions that 
if the adiabatic exponent 7 is larger than 3/2 and there exists a unique steady state, then the 
density of any global weak solution converges to the steady state density in some L 9 -norm as 
time goes to infinity. Under the assumptions that the adiabatic exponent 7 is close to 1 and 
the external forces satisfy some decay properties in the far field, Matsumura and Yamagata [27] 
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studied the global existence and large-time behavior of weak solutions to the Cauchy problem of 
(jl.ip . (|1.2p when the initial perturbations are suitably small in 1? n L°° for density (away from 
vacuum) and in H 1 for velocity. Recently, under the same smallness assumptions on the initial 
perturbations, Li and Matsumura [22] succeeded in removing the smallness condition on | -y — 1 1 
and the decay assumptions on the external force, and thus, improved the Matsumura- Yamagata's 
result [27]. 

It is worth noting that among the papers |10 p i H l22 j l27| mentioned above, the initial pertur- 
bation of density around a given positive state is small in L°°, which in particular implies that 
the density is uniformly away from vacuum. However, as emphasized in many papers related 
to compressible fluid dynamics [2TH1I35]. the possible presence of vacuum is one of the major 
difficulties when the problems of global existence, uniqueness and regularity of solutions to the 
compressible Navier-Stokes equations are concerned. 

Thus, our main aim in this paper is to establish global well-posedness theorem and to study 
the large-time behavior of solutions for the compressible Navie-Stokes equations (jl.ip . (11 .2[> with 
large external forces and vacuum states. 

Before stating the main results, we explain the notations and conventions used throughout 
this paper. For simplicity, we denote 

fdx = / fdx. 

For 1 < p < oo and integer k > 0, we adopt the following simplified notations for the standard 
homogeneous and inhomogeneous Sobolev spaces: 



LP = LP(R 3 ), W k ' p = W k ' p (R 3 ), H k = W k ' 2 , 
D 1 = {«£ L 6 (M 3 ) | ||Vu|| L 2 < oo} , 
k D l,p ={«£ Aoc( m3 ) I I|Vw||xp < oo} . 

We first study the stationary problem of (jl.lj) - (jl.5j) . In view of |22[ Remark 2.1], it is 
known that the smooth steady solution, (p s (x), u s (x)), is unique and u s = 0. Hence, we infer 
from (|1.2|) that the steady state density p s = p s (x) satisfies 

VP(ps) = p s V/, p s {x) ->■ poo as \x\ ->• oo, (1.6) 

which implies that p s is uniquely determined by 

" S p- 1 P'(p)dp = f{x). 



Po 



In order to avoid the vacuum states in p s , we suppose that 

fPoo pii p \ fOO pit N 

- / -^-dp < inf f(x) < sup f(x) < / -^dp. (1.7) 

Jo P zelR 3 xgM 3 ho* P 

Thus, it follows from (fL6j) and ([L7]) that 

Proposition 1.1 Assume that f G H 2 satisfy {1. 7\ ). Then the stationary problem il.b]) has a 
unique solution p s = p s {x) satisfying 

Ps~ Poo G H 2 , < p < inf p s (x) < sup p s (x) < p < oo, (1.8) 
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where p,p are two positive constants depending only on A, 7,Poo, inf f(x), and sup fix). Fur- 
thermore, if f E W 2,q with some q E (3,6), then 

l|Vp s || < C, (1.9) 

where C is a positive constant depending only on A, 7, p^, inf f(x), and \\f\\ H 2 nW 2, g . 

x&R 3 

Remark 1.1 For P(p) = Ap 1 with A > 0,7 > 1, it is easily seen from hi. 6\) and fli. 7| ) that 

1 

7-1 , 7-1 

Hsy-^j — \ p 

t2 



ps{x) = up- 1 + -^-m 



which, together with f E H , implies that il.8\) holds provided 

inf f( x)> -^Lpl^ 



7-1 
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Remark 1.2 To study the large-time behavior of global weak solutions, the authors in \22\ 21 ^ 
technically required that f E which particularly yields Vp s E L°° and plays a key role in their 
analysis. Here we only assume that f E H 2 which is much weaker than that in \22l\21^ and is 
only used to guarantee hi. 8) . The additional condition f E W 2,q D H 3 will be used to derive the 
high order estimates needed for the global classical existence. 

With the steady state (p s ,u s ) at hand, we can define the initial energy Co as follows: 

C = J^(G(p ) + ^ Po \u \ 2 ^dx, (1.10) 
where G(-) is the potential energy density defined by 

J p s J p s S J p a S 

Now, our first result concerning the global existence of classical solution of (|l.ip - (|1.5|) can be 
formulated as follows. 

Theorem 1.1 For q E (3, 6), let f E H 2 n W 2 ' q satisfy JJ.Tp and p s = p s (x) be the steady state 
density of hl.b}) . For given positive numbers M (M may be arbitrarily large) and p > p + 1, 
assume that the initial data {pq,uq) satisfy 

(po-Poc,P(po)-P(Poc))£H 2 nW 2 > q , u £H 2 , (1.12) 

< inf p (x) < sup p (x) < p, ||Vu |||a < M, (1.13) 

and that the compatibility condition 

- pAu - (A + /x)Vdivu + VP( P o) = pl /2 g (1.14) 

holds for some g E L 2 . Then there exists a positive constant e > 0, depending only on p, X, A, 
7; Poo, P, M, inf f(x) and ||/|| H 2 r\W 2 >i> such that if the initial energy satisfies 

Co < e, (1.15) 
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the Cauchy problem A1.1\ )- ITL5\) has a unique global classical solution (p, u), defined on M 3 x (0, T] 
for any < T < oo and satisfying 

0<p(x,t)<2p for all x G R 3 , t > 0, (1.16) 

and 

' (p - Poo ,P(p) - P(poo)) G C([0,T]; H 2 n W 2 ><0, 

< u G C([0, T]; -ff 2 ) Pi L°°(t, T; H 3 n l^ 3 ' 9 ), (1.17) 
k «t G L°°(r , T; U 2 ) fl H^(t , T; H 1 ), 
for any < r < T < oo. Moreover, one has the following large-time behavior: 

lim (||p(-,i) - Ps \\ LP + ||Vu(-,i)|U0 = 0, (1.18) 

t— >-oo 

TOt/i any 

p€(2,oo), re [2,6). (1.19) 

Remark 1.3 It is clear from ji. 17\ ) that the solution obtained in Theorem \l.l\ becomes a classical 
one for any positive time ( 11 5[ Remark 1.1]). Moreover, although the solution has small energy, 
its oscillations could be arbitrarily large and the interior vacuum states are allowed. 

Remark 1.4 Recently, in the absence of external force, Huang-Li-Xin ]i7f (see also ^]) proved 
the global existence of classical solutions of the Cauchy problem of \1.1\) . hi. 2) with smooth 
non-negative initial density under the following compatibility conditions 

- pAu - (A + /i)Vdivu + VP(po) = Po9 (1-20) 

with g = g(x) satisfying 

pl /2 g G L 2 , Vg G L 2 , (1.21) 

which play a key role in the analysis of f4\\ 17$. It is worthwhile noting that the compatibility 
conditions kl.20\) . $1.21)) are much stronger than that in (O^ ). 

It is well-known that the discontinuous solutions (namely, weak solutions) are fundamental 
and important in both the physical and mathematical theory. So, our second aim is to study 
the global weak solutions (see Definition 1 1 . 1 1) of (|l.ip - (jl.5p . 

Definition 1.1 A pair of functions (p,u) is called a weak solution of li 1. 1\) - [T75\) , provided that 

P-Poo G^c(0,oo;L 2 nL°°), u G L^ c (0, oo; H 1 ), 
and that for all test functions ip G D(IR 3 x (— oo, oo)) and j = 1, 2, 3, 

//"OO /* 
po f ijj(^0)dx + / (pipt + pu • V^O dxdt = 0, 

p uiip(-,0)dx + J J (pu j ^ t + pu j u ■ Vtp + P{p)il) Xj ) dxdt 

poo r 

(pVv? ■ Vip + (p + X)(divu)^ Xj ) dxdt — / / pf X] ipdxdt. 
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The global existence of weak solutions of (ll.ip -( fl~5j) with discontinuous initial data can be 
stated as follows. 



Theorem 1.2 Let f G H 2 satisfy fli. 7| ) and p s = p s (x) be the steady state density of lll.6\) . For 

given positive numbers M (M may be arbitrarily large) and p > p + 1, assume that the initial 
data (pq,uq) satisfy 



(p ~ Poo, P{po) ~ P(Poo)) G L 2 n L°°, n G H 1 , 
< inf po(x) < sup po(aO < Pi ||Vtto|| r 2 < 

^GM 3 xeR 3 



Cl-22) 



T/ien i/iere exists a positive constants > 0, depending only on p, X, A,j, p^, p, M, inf /(x), and 

x£R 3 

||/||^2, such that if 

C <e, (1.23) 

the Cauchy problem \1. 1\) - [1~5\) has a global weak solution (p, u) on M 3 x (0, oo) in the sense of 
Definition li.il satisfying 

0<p{x,t)<2p for all x G R 3 , t > 0, (1.24) 
p-/>oo G C([0,oo);L 2 ), pn G C([0, oo); F" 1 ), Vn G L 2 (0, oo; L 2 ) (1.25) 
and /or any p G (2,oo), 

lim (||p(-,t)-p s |U P + ||'u(-,t)|U J .ni-) = 0. (1.26) 

t— >oo 



Remark 1.5 Theorem li.j?l extends those results in J1(M1^[2R 271 to the case that both the 
vacuum states and the large external forces are involved. Moreover, the regularity condition 
f G H 2 on the external force is much weaker than the one f G H 3 which is technically needed 
in I22[ \27y. Indeed, more regularities of the solutions away from t = can be obtained (cf. \lCfj )- 

Remark 1.6 Similar to [11], the condition ||Vno|| 2 2 < M in both \1.13\) and $1.22\) can be 

replaced by Hnoll^o < M with any j3 G (1/2, 1]. 

We now comment on the analysis of this paper. Note that the compatibility condition 
(|1.14p is much weaker than those in [3], i.e., (|1.2U|) . (|1.21|) (see Remark ll.4|) . we cannot apply 
the local existence theorem of classical solution in [4] to the problem considered. Indeed, we shall 
split the proof of Theorem 11.11 into three steps. Roughly speaking, we first use the well-known 
Matsumura-Nishida's theorem (see Lemma 12. 5p to guarantee the local existence of classical 
solutions with strictly positive initial density, then extend the local classical solutions globally 
in time just under the condition that the initial energy is suitably small (see Proposition I5.ip . 
and finally let the lower bound of the initial density go to zero. So, to this end, we need some 
global a priori estimates which are independent of the lower bound of density. It turns out that 
the key issue in the proof is to derive both the time-independent upper bound of density and 
the time-dependent higher order estimates of (p,u). To do this, we will borrow some ideas due 
from |10 p i H [T7 1 l22j. However, because of the arbitrariness of external forces, the presence of 
vacuum states and the weaker compatibility condition (|1.14p . some new difficulties arise and the 
methods therein cannot be applied directly. 
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First, similar to that in (HUE], we begin our proof with the careful initial layer analysis. 
To do this, we technically need the following modified "effective viscous flux": 



F 4 pji [(A + 2p)dwu - (P(p) - P{p s ))] , 



(1.27) 



which was introduced by Li-Matsumura [22] and is different from the ones in |10 p il | [r7 ] l27|. Bas- 
ing on (|1.2|) . (jl.27p and the standard L p -estimate of elliptic system, one can derive some subtle 
connections among the modified "effective viscous flux", the gradient and material derivative 
of the velocity, and the pressure (see Lemma l3.3p . which are important in the entire analysis, 
particularly in closing the time-independent energy estimates stated in Lemma 13.21 To obtain 
such connections, we need to deal with some difficulties induced by the large external forces. 
This will be done by making a full use of the mathematical structure of the steady state density 
(see (ll.6j) ) and adopting an idea due to Huang-Li- Xin [16] (see (|3.29p below), which enable us 
to control the terms associated with the pressure and the large external force by the deviation 
of the density p from the steady state p s . Moreover, compared with [TUILT7] , the properties of 
||V-F||£2 and ||Vcurlti||^2 are proved in a different manner due to the fact that f,p s only belong 
to H 2 , which also makes the analysis here need to be more careful than that in |22j . 

However, unlike that in [22], it seems difficult to use directly this modified "effective viscous 
flux" F in (jl.27p to prove the uniform upper bound of the density p since we only assume that 
/ G H 2 which implies that Vp s G L s with any s G [2, 6]. Indeed, we overcome this difficulty by 
replacing F by the following standard "effective viscous flux" (see Lemma |3.8() : 



which is in a similar form as the one defined in |10 p il t fT7 l l24j. A key observation is that for 
r, v\ G (3, 6], T2 G (6, oo) and 1/rx + \jri = 1/r, one has 



where " ' " denotes the material derivative. Estimate (jl.29p shows that we only need that p s 
satisfies Vp s G L 2 n L 6 . This is the main difference between F and F defined respectively in 
(jl.27p and (jl.28p . By virtue of (jl.29p . we can apply the Sobolev embedding inequality to derive 
a desired estimate of ||-F||l°°) an d thus, prove the pointwise boundedness of density by using the 
Zlotnik inequlaity (see Lemma l2.2p . 

There are also some new difficulties lying in the proof of the higher order time-dependent 
estimates. Indeed, to achieve the estimates on the derivatives of the solutions, we first prove the 
important estimates on the gradients of the density and velocity by solving a logarithm Gronwall 
inequality in a similar manner as that in [171118] , As a result, one can also easily obtain the 
L 2 -estimates for the second-order derivatives of density, pressure and velocity. However, due to 
the weaker compatibility condition in (11 . 14H (cf. (I1.20|) . (11.2ip ). the method used in [17] cannot 
be applied any more to obtain further estimates needed for the existence of classical solutions. 
In fact, instead of the L 2 -method, we succeed in obtaining these classical estimates by deriving 
some desired L 9 -estimates (3 < q < 6) on the higher-order time-spatial derivatives of the density 
and velocity, basing on some careful initial-layer analysis (see Lemmas 14.41 - 14,61) . 

The rest of this paper is organized as follows. In Sect. 2, we first recall some known 
inequalities and facts which will be frequently used in our analysis. In Sect. 3, we derive the key 
a priori estimates of the weighted estimates on the gradient and the material derivative of the 



F 4 (A + 2p)divu - (P(p) - P( Ps )) 



(1.28) 



VF\\ L r < C {\\ P u\\ L r + \\{ P - Ps )Vf\\ L r) 

KCiWpuhr + WVfhnWip-p^Ur,) 



(1.29) 
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velocity and the pointwise upper bound of the density, all of which are independent of t and will 
be used to proved the large-time behavior. In Sect. 4, we prove the time-dependent estimates on 
the higher-order norms of the solutions, which are needed for the existence of classical solutions. 
Finally, the main results (i.e. Theorems 1 1 . 1 1 and 1 1 . 2 p will be proved in Sect. 5. 

2 Auxiliary lemmas 

In this section, we list some elementary inequalities and known facts which will be used 
frequently later. We begin with the well-known Gagliardo-Nirenberg-Sobolev-type inequality 
(see [20]). 

Lemma 2.1 Let f G D 1 and g G L q nD 1,r with q G (l,oo) and r £ (3,oo). Then there exists a 
positive constant C , which may depend on q and r, such that 

WSWifi < ciiv/iu (2.i) 

\\9\\V < C\\ 9 ||^-3)/(3r+ 9 (-3)) n Vg ||3 : /(3r+ g (r-3))_ ^ 

The following Zlotnik inequality, whose proof can be found in [36], will be used later to 
prove the uniform-in-time upper bound of the density. 

Lemma 2.2 Assume that the function y G W 1,l (Q,T) solves the ODE system: 

y' = g(y)+b'(t) on [0,T], y(0) = y , 

where b G W^^T) and g G C(R). If g{oo) = -oo and 

b(t 2 ) ~ Kh) < N + JVi(t 2 - h) (2.3) 

for all < t\ < t2 < T with some positive constants N$ and N\, then one has 

y(t) < max{y ,O + < +oo on [0,T], (2.4) 

where ^* G R is a constant such that 

9(0<-Ni for e>T- (2-5) 

In order to obtain the time-dependent estimates of ||Vu||£°° and ||Vp||i2 ni 6, we need the 
following Beale-Kato-Majda-type inequality, the proof of which can be found in [18J. 

Lemma 2.3 For 3 < q < oo, assume that Vti 6l 2 H D 1,q . Then there exists a constant C > 0, 
depending only on q, such that 

||Vu||loo < C (||divu||i<x> + ||V x u||loo ) In (e + || V 2 u\\ Lq ) + C||V«|| L 2 + C. (2.6) 

Since there is no vacuum state in the far field, it is easy to show that u G L 2 even that the 
density is only nonnegative. Indeed, we have 

Lemma 2.4 Let p s = p s (x) be the steady state density as in Proposition \l.l\ Assume that (g,v) 
satisfies the following conditions: 

0<g<2p, (q- p s ,0 1/2 v) G L 2 , VveL 2 . 

Then there exists a positive constant C, depending only on p, p and p, such that 

\Mv < C (\\Q 1/2 v\\ L 2 + \\q- psf^WVvW^ . (2.7) 
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Proof. Indeed, it is easy to see that 



J \v\ 2 dx<C(p) J p 2 \v\ 2 dx <C(p) (J g 2 \v\ 2 dx + J (g - p s ) 2 \v\ 2 dx 
<C(p,p) j g\v\ 2 dx + C(p) (J \g-p s \ 3 dx^j ' (J \v\ 6 dx^j 



< C{pr P )\\g l/2 v\\h+ C(p,p,p)\\g-p s \\% 3 \\Vvf L2 , 

which immediately proves (|2.7h , □ 
We end this section with the following local well-posedness theorem of classical solution to 
the problem (ll.ip - (11.5j) when the initial density is strictly away from vacuum (see, e.g. [29], and 
especially Matsumura-Nishida \25\ Theorem 5.2]). 

Lemma 2.5 Assume that the initial data (po,uo) satisfies 

(po ~ Poo, «o) e H 3 , inf p Q (x) > 0. (2.8) 

Then there exist a small time Tq > and a unique classical solution (p, u) to the Cauchy problem 
COP-COP on M. 3 x (0, T ] such that 

inf p(x,t) > - inf p (x), (2.9) 

and 

(p-PooeCdO^H^nCH^To^H 2 ), 

[ueCi^To^H^ncH^To^H^n^^To-H^, <2 ' l0) 

where Tq > may depend on inf po(x). 



3 Time-independent a priori estimates 

This section is concerned with the time-independent (weighted) energy estimates and the 
uniform upper bound of density, which are essential for the proofs of Theorems 11.11 and ll.2i To 
do this, we assume that (p, u), defined over (0, T) with some positive T > 0, is a smooth solution 
of the Cauchy problem (jl.ip ~ (jl.5p . For simplicity, we introduce the following functionals: 

$i(T) = sup a [ \Vu\ 2 dx+ [ a [ p\u\ 2 dxdt, (3.1) 

0<t<T J JO J 

$ 2 (T) = sup a 3 [ p\u\ 2 dx+ [ a 3 [ \Vu\ 2 dxdt (3.2) 

0<t<T J Jo J 

and 

' |2 



$ 3 (T) = sup / \Vu\ 2 dx, (3.3) 

0<t<T J 

where a(t) = {1, t}, and the symbol " ' " denotes the material derivative v = vt + u ■ Vv. 
The main purpose of this section is to prove the following key a priori estimates. 
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Proposition 3.1 Assume that the conditions of Theorem \1.S\ are satisfied. There exist two 
positive constants e and K, depending only on u, A, A, 7, p, p, p, inf fix), WfWu 2 an d M, such 

that if (p, u) is a smooth solution of hi. 1\) - [T75\) satisfying 



[0 < p(x,t) < 2p for all (x,t) el 3 x [0,T], 
1 $i(T) + $ 2 (T) < 2C 1/2 and ^ 3 (a(T)) < 3K, 

then one has 

( < p(x, t) < \p for all (x, t) G M 3 x [0, T] , 
1 $i(T) + $ 2 (T) < C 1/2 and * 3 (*(T)) < 2K, 
provided that the initial energy Cq defined in il.lO\) satisfies 



(3.4) 



(3.5) 



C < e. (3.6) 

Proof. Proposition 13.11 follows directly from the following Lemmas 13.51 13.61 and 13.81 with 
K and e being the positive constants as in Lemmas 13.51 and 13. 8| respectively. □ 

Remark 3.1 We assume throughout this section that the initial data and the external force only 
satisfy the conditions of Theorem \1. 6 A and hence, the uniform-in-time estimates derived in this 
section can be used to study the existence and large-time behavior of global weak solutions as 
stated in TheoremM.SX 



For notational convenience, throughout this section we denote by C or Cj (i = 1,2,...) the 
generic positive constants which may depend on p,,\,A,j,p,p,p, inf f(x), WfWn 2 an d M, but 

not on T. We also sometimes write C(a) to emphasize that C relies on a. 
We start the proof with the following standard energy estimate. 



sup / (-p\u\ 2 + G(p)]dx+ [ [ (p\Vu\ 2 + (p + X)(divu) 2 ) dxdt < C , 

t<t<TJ V 2 J Jo J 



(3.7) 



iivn, 



Lemma 3.1 Let (p,u) be a smooth solution to il.l\) - [T75\) on M 3 x (0, T]. Then, 
I [t-p\ u \ 2 + G{p))dx+ f 1 

0<t<T 

where G(p) is the potential energy density defined in M.ll}) . 
Proof. Thanks to (|1.6|) . the momentum equation (|1.2|) can be written as 

pu t + pu-Vu+ (VP{p) - pp~ l VP{p s )) = pAu + (p + A)Vd 
which, multiplied by u and integrated by parts over M 3 x (0,i), yields 

l - J p\u\ 2 dx + J u- {VP(p) - pp^VPips)) dxds 

+ J J (p\Vu\ 2 + (p + \)(diYu) 2 )dxds = 0. (3.8) 
After integrating by parts, one infers from (jl.ip that 

(VP(p) - pp^VPips)) dxds 
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- J q J div(pu) (^j P ^-jr-<%) dxds 

Jo I pt {J P ^T~ d ^) dxds = I Gi - p)dx 







which, inserted into (|3.8p . leads to the desired estimate in (|3.7p . □ 
It is clear that for all < p < 2p and p < p s < p, there are positive constants C\ and C2 
depending only on p, p, and p, such that 

Ci(p-p s ) 2 <G(p)<C 2 (p-p s ) 2 , 
so that, it readily follows from (|3.7p that 

sup \\p-psWh < CC . (3.9) 

0<t<T 

The next lemma is concerned with the temporary weighted L 2 -estimates on the gradient 
and the material derivatives of the velocity, the proof of which will be concluded in Lemma 13.61 
below. The idea of the proof mainly comes from |1Q |, I1T ^ [22] . However, due to the arbitrariness 
of external potential forces, the presence of vacuum states, and especially, the weaker regularity 
assumption of the external forces (i.e. / 6 H 2 ), the analysis here needs to be more careful. 

Lemma 3.2 Let (p,u) with p £ [0, 2p] be a smooth solution of lil. 1) )- [L5\) on IR 3 x (0, T]. Then 
there exists a positive constant C , depending on p, such that 

$i(T) < C ( C o + ^ a\\Vu\\l 3 dt^ , (3.10) 

M T ) < C (c + $i(T) + a 3 \\Vu\\l 4 d?j . (3.11) 



Proof. It follows from ([L2]) and {T5]) that 

pii - pAu -(p + A)Vdivn + V (P(p) - P(p s )) = (p - p 8 )Vf. (3.12) 
For m > 0, multiplying (|3.12p by a m u in 1? gives 

J p\ii\ 2 dx = pa m j Au ■ iidx + (p + A)a m J Vdivu • iidx 

-a m Ju-V (P(p) - P{p.)) dx + a m J(p- p s )u ■ Vfdx 
4 



i=l 



where the right-hand side can be estimated term by term as follows. First, by the definition of 
the material derivative " ' " and integration by parts, we easily get that 

Ji = -| {a m \\Vu\\l 2 ) t + ^ma^a'WVuf^ - pa m J dm 3 di(u k d k u j )dx 
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_p 

and similarly. 



< ~ (^ m ||Vn||| 2 ) t + Cma m - x a'\\Vu\\l2 + Ca m \\Vu\\ 3 L3 , (3.14) 



h < "^y^ {° m \\&M\h) t + CW m -V||Vu||| 2 + C(7 m ||Vn|||3. (3.15) 
Secondly, noticing that (jl.ip implies 

P(p)t + div(P(p)u) + (pP'(p) - P(p))divu = 0, (3.16) 
so that, using fjl .8[) . (|2.ip and (|3.9p . we obtain after integrating by parts that 



/ 3 = / k m (divn) 4 (P(p) - P(p s )) - a m (u • Vu) • V (P(p) - P(p s ))] dx 

= j t J <r"Mivu (P(p) - P{ Ps )) dx - ma m - l a' J divu (P(p) - P(p s )) 

+ y cr m ((pP'(p) - P(p))(divu) 2 + P(p)d i u j d j u i + u ■ Vu ■ VP(p s )) dx 

<jj a m divu(P(p)-P(p s ))dx + ma m - 1 a'\\P(p)-P(p s )\\ L 4Vu\\ L2 
+Ca m \\Vu\\\ 2 + Ca m \\u\\ L 4\Ju\\ L 2\\V Ps \\ Li 



d 

< — 
~ dt 



J a m dwu{P{p) - P{ Ps ))dx + C|| VtiUla + CroV^-VCb. (3.17) 



Analogously, using the fact that pt + div((p — p s )u) + div(p s u) = due to (jl.ip . we have from 
Ob. dO and (ESI) that 



/ 4 = A J a m( p _ ps } u . V f dx _ ma m-l a l J^ p _ ps)u . yj dx 

+a m J [div(p s u)u ■ V/ - (p - Ps )u ■ V(u • V/) + (p - p s )u • Vu ■ V f] dx 
-Jtl ° m (P~ Ps^-Vfdx + Cma^a'llp- p s \\ L 2\\u\\ L 4Vf\\ L s 
+Ca m J (|n||Vn||V/| + |Vp s ||n| 2 |V/| + |p - p s ||u| 2 |V 2 /|) dx 

-Jtl - • v / da; + c^ m - Vc^^iiv^H^ + C^H^lUa ||V^|U 2 ||V/|| i3 

+Ca m ||Vp s || L 3||n|| 2 , 6 ||V/|| L 3+Ca m ||p-p s || i6 ||n|| 2 :6 ||V 2 /|| L2 

-Jtf ~ Ps)M ' Vfdx + C||Vu| ^ 2 + C™ 2 ° a{m ~ 1) ° , CQ' ( 3 - 18 ) 
Combining (^T3]) - (l3TT5l ). (I3TT7D and (I3TT8D leads to 

^ m ||Vn||l 2 + ^±^a m ||divn||l 2 ) i + J a™p\u\ 2 dx 

J a m (divn (P(p) - P(p s )) + (p - p s )u • V/) 
+CmV( m ~ 1 VQ ) + C(l + mcr^V) || VuH^ + Ca m ||V-u||| 3 , (3.19) 



d 

< 
~ dt 
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where the first term on the right-hand side can be estimated as follows: 



a m (divu(P(/>) - P(p s )) + (p- Ps )u ■ Vf) dx 



< Ca m (Hp - /9 S || L 2||v«|| L 2 + ||p - Ps \\lAHl« II v/|| l3 ; 

< ^a m \\Vu\\ 2 L2 + Ca m C . 



Thus, choosing m = 1 in (|3.19p . integrating it over (0, T) and using (|3.Tj) one gets (|3.10p . 

To prove (|3.1ip . operating a m u 3 (dt + div(u-)) to both sides of (|3.12p and integrating the 
resulting equations over IR 3 , we obtain after summing them up that 



i=i 



d t (pu j ) + d k (pu k u j )] dx - pa m \ v? \a U { + d k (u k Au j 



dx 



(p + A)cr m / ii J djdhrut + d k (u k djdivu) 



dx 



\djP{p)) t + d k (u k d, (P( P ) - p(ps)))] dx 

2 

Pt djf + d k (u k (p- p s )djfj dx = ^Ri. 



(3.20) 



i=i 



We now estimate each term in (|3.20p . First, by (II. lh one easily gets that 

1 d 



Li 



r> ^ J a m p\u\ 2 dx - ^<r m - V J p\ii\ 2 dx. (3.21) 
Recalling the definition of " ' " , we deduce from integration by parts that 
L 2 = pa m J (\Vu\ 2 - d k ii j d k (u ■ Vu j ) + d k u j u k Au j ^j dx 

= pa m J (\Vu\ 2 - d k u j d k u l d lU j - d k v?u l d 2 kl v? - d 2 kl v?u k d { v? - d k v? d x u k fyu^ dx 

- ^ /(iv*|»- o t M„ + 3M,y - at ,W) * 

\Vu\ 2 dx - Co m J \Vu\ A dx, (3.22) 
where the Cauchy-Schwarz inequality was also used in the last inequality. Analogously, 
L 3 > (p + A)cr m J (diviifdxds - C o m J \Vu\ \Vu\ 2 dx 



Gu + A)(divn) 2 -^\Vu\ 2 )dr -Crr' 



In view of (|3.16p . we obtain after integrating by parts that 

Ri = -o m J divu (div(P(p)w) + (pP'(p) - P{p))d 



\Vu\ dx. 



ivu) dx 



(3.23) 



o m \ [u k d k d\vu + d k u j d jU k ) (P(p) - P(p s )) dx 
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= -a m j divu [(pP'(p) - P(p))divu + div (P(p s )u)] dx 
-a m J d k v?d jU k (Pip) - P(p s )) dx 

< g* m [|V«[|£i + Ca m {\\Vuf L2 + \\uf L6 1| Vp s f L s) 

< ^ m ||Vu||| 2 + Ca m \\Vu\\ 2 L2 , (3.24) 
8 

where we have also used (|1.8p and ()2.ip . Similarly, by (II. ip we have 

R 2 = a m J pu-V (u j djf) dx-a m j d k ii j (u k (p - p s )d j f^j dx 

< Ca m (||v«|| i2 ||v/|U 3 |HU6 + HHi^iMMivVlM 

< |a m ||Vn||| 2 + Ca m (\\p l/2 u\\ 2 L2 + ||Vu||^) , (3.25) 
where we have also used Cauchy-Schwarz inequality and the following estimate: 

II -II ^ 1/2 ■ 1 1 3-/2 1 1 ■ II 1/2 ^ n \\ 1/2 - iil/2||„ . ,,1/2 

\\pu\\ L 3 < C\\p ' u\\£ 2 \\u\\£ 6 < C\\p 1 u\\£ 2 \\Vu\\£ 2 . 
Putting the estimates of (ET2"T]) - (pI^5l into (I3T2T1 gives 

[a m \\p l ' 2 u\\l 2 ) t + pa m \\Vu\\l 2 
< Ca rn (||V«||| 2 + ||Vu||£ 4 ) + C (ma^a' + a m ) \\p 1/2 u\\ 2 L2 . (3.26) 

Thus, choosing m = 3 in (|3.26p and integrating it over (0, T) lead to (|3.1ip immediately. □ 
Next, we prove some important connections among the modified "effective viscous flux", the 
density, and the gradient and material derivative of the velocity, which are crucial for our further 
analysis. To do this, we need to make a full use of the mathematical structure of steady state 
density p s of (|1.6p and the standard L p -estimate of elliptic system, basing on a key observation 
due to Huang-Li-Xin [16] (see (|3.29j) below). 

Lemma 3.3 Let f £ H 2 satisfy JJ.7| ) and p s = p s (x) be the steady state density of hl.b}) . 
Assume that (p, u) with p G [0, 2p] is a smooth solution of \l. 1\) - (T75\) and F = F(x, t) is the 
modified "effective viscous flux" defined by {1.21 ). Then there exists a positive constant C, 
depending only on p,p,p, inf f(x) and ||/||#2, such that 

||VF|| L2 + HVO^curlu)!!^ < C (\\pu\\ L 2 + ||V«|| L a + \\p - p s \\ 2 L6 ) , (3.27) 
||Vti|| L 6 < C (\\pu\\ L 2 + ||Vn|| i2 + \\p - p s \\ 2 L e + \\p - Pshe) ■ (3.28) 

Proof. The proofs are analogous to those in |10pi7j. however, since p s — Poo H 3 , we modify 
the proof of ||VF|| I/2 slightly. Indeed, as observed in [16J, one can utilize (jl.6p to get that 

p- 1 (VP(p) - P Vf) 

= pj 1 [V (P^) - P( Ps )) - p- l (p - P s)VP(Ps)] 

= V [p; 1 (P(p) - P(p s ))] + p- 2 [Pip) - Pip,) -ip- Ps )P'iPs)] V Ps . (3.29) 
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Thus, by virtue of ([T5j> . (fL2Tl) and ([229]), we deduce 



A 



which gives 
and 



p s pii — VF + /icurl (p s 1 cur hi) 

= - [(A + 2/x)(divu)Vp7 1 - pVpJ 1 x (curlu)] 
+ [P{p) - P{p s ) - P'{ Ps )(p - p s )] Vp- 1 

G 1 +G 2 , (3.30) 
AF = div(/57 V -Gi- G 2 ) (3.31) 



pA (p s curlu) = /iVdiv (p s 1 curln) — /icurlcurl (p s 1 curln) 

= pV (curlu • V/^ 1 ) - curl^/w - G x - G 2 ). (3.32) 



Noticing that 



IIGilUa + HGalUa <C(||Vp s || L 6||Vn|| L 3 + ||V /3s || L 6||( /9 -p s ) 2 || i 3) 
<C(\\Vu\\ L 3 + \\p-p s \\ 2 L6 ), 

we apply the standard L 2 -estimate to (|3.3ip and (|3.32p to obtain 

llVFH^ + HVO^curWHI^ 

< C(\\pu\\ L 2 + \\Gx\\ L 2 + \\G 2 \\ L 2 + ||VpJ r 6||Vu|Us) 

< C (\\pu\\ L 2 + \\Vu\\ L 3 + \\p - PsWle) , 

which gives (|3.27p . 

On the other hand, using (fL27jl . (l2TTjl . and (137271) . we find 

HVnll^e < C (\\divu\\ L 6 + HcurlnH^) 

< G (\\F\\ L 6 + \\p - p s \\ L 6 + Hp^curl-ull^e) 

< C {\\VF\\ L 2 + \\p - p s \\ L e + HVC^cuilu)!!^) 

< C(||H| L 2 + llp-Psllle + \\p~ Ps\\l* + l|Vu|| L 2) + -||Vn|| L 6, 

where we have used the fact that ||Vu|| 2 3 < C\\ Vn||^2 ||Vu||^6 and Cauchy-Schwarz inequality 
in the last inequality. This proves (13.28|) . and thus, the proof of Lemma 13.31 is completed. □ 
With the help of (j3.4|) and Lemma 13.31 we can prove 

Lemma 3.4 Let (p,u) with p G [0, 2p] be a smooth solution of n.l\) - fL5\) on R 3 x (0, T]. Then 
there exists a positive constant Eq > 0, depending on p, such that 

[ a 3 (\\Vuf L4 + Wp-psWi* + \\F\\ L4 + ||curht||£ 4 ) dt < CC . (3.33) 
Jo 

provided $i(T) + $ 2 (T) < 2C 1/2 and C Q < e . 
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Proof. Direct calculations from (jl.ip and (ll.27|) show 



(p " Ps)t+ vh (p(/,)_py) 



-u ■ V(p - p s ) - (p- p s )divn -u-Vp, 



s 2p + X' 

Multiplying this by 4(p — p s ) 3 in L 2 and integrating by parts, we find 

^ J( P - Ps ) A dx + — ^ J Ps (P( P ) - P( Ps )) (p - p s fdx 

<C f {{ P -p s f\Vu\ + \p-p a \M\Vps\ + \p-ps?\F\)dx 

< C\\p - p s \\l, (||Vn|| i2 + ||n|| i6 ||Vp s || L3 ) + C\\p - p s \\U\F\\^ 
<6\\p-p s \\i, + C(6)(\\Vu\\l 2 + \\Ff Li ), 5>0, (3.34) 

where we have used (jl.8p . (|2.ip and Cauchy-Schwarz inequality. Due to (jl.8p . one has 

p s (P(p) - P( Ps )) (p - Ps f = (p- p s fp s f P'(ap + (1 - a) Ps )da 

Jo 

>C{ P - Ps f 

for some positive constant C depending only on A, 7, p, p and p. Hence, choosing 5 > suitably 
small, multiplying (|3.34p by <r 3 and integrating it over (0,T), we infer from (|3.7p and (|3.9p that 

^^Hp-p.H^^^CCo + C [ T a 3 \\F\\ 4 L4 dt. (3.35) 
Jo Jo 

By (fL27]) . we have 

||Vn|| i4 < C(||div U || L 4 + ||curlw|U*) < C (\\F\\ L * + \\p - p s \\ L , + ||p7 1 curl«|| L 4) , 
which, combining with ()3.35|) and (|3.27|) . yields that 

l-T 



[ a 3 (||Vn||| 4 + ||p-p s ||i 4 + ||F||l 4 )^ 
Jo 

<CC + C [ a 3 (||F||| 4 + H^curlulli*) dt 
Jo 



T 



<CC + C ^(llFll^llVFlll. + ll^curlnll^llV^curMH 3 ,)^ 
Jo 

<CC + C^ a 3 (\\Vu\\ L 2 + C 1/2 ) (|H|| 2 + ||Vu||£a + \\p - Ps \\%) dt 

< CC + c£ (o- 1/2 \\Vu\\ L2 +C] /2 ) {a 3 \\pu\\ 2 L2 ) 1/2 a\\puf L2 dt 

+ c£ \\Vu\\ L2 (JvWWVuWv) ^\\Vu\\ 2 Li + a 3 Cl /2 \\p - p s \\%) dt 

+ CC]' 2 f T a 3 (||Vu||| 2 + ||Vn||| 4 + \\p - p s \\%) dt 
Jo 

< cc + cc]' 2 [ a 3 (||Vu||£ 4 + ||p - p s \\%) dt, 
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where in the last inequality we have used $i(T) + $ 2 (T) < 2C / . This directly gives (|3,33p 
provided C < e = min{ 1 , (2C) " 2 } . □ 

Lemma 3.5 Let (p,u) with p 6 [0, 2p] 6e a smooth solution of lll.l\ )- {T75]) on M 3 x (0, T]. Then 
there exist positive constants K and E\, depending on p and M , such that 

ra(T) 

$ 3 (a(T)) + / \\p 1/2 u\\ 2 L2 dt<2K, (3.36) 
J o 

provided $ 3 ((r(T)) < 3K and C < £i- 



Proof. Choosing m = in (|3. 19[) and integrating it over [0,<r(T)], we deduce from (|3.7p . (|3.9 
(I3.28P and Cauchy-Schwarz inequality that 



/•°v ) 

$ 3 (a(T)) + 2 / \\ p y 2 uf L2 dt 
Jo 

,<x(T) 

<C(C + M) + C / ||V-u|| 3 3 dt 
Jo 

< C (Co + M) + C / HVnll^llVnll^ 2 ^ 
Jo 

f a(T) ra{T) 

<C(C + M) + C \\Vu\\ 6 L2 dt+ Wp^iif^dt 
Jo Jo 

ra(T) 

<K + CC [<*> 3 (a(T))] 2 + / \\p 1/2 u\\ 2 L2 dt 

Jo 



with a positive constant K = C{Cq+M) depending on p, Cq and M. As a result, we immediately 
obtain (1336]) provided $ 3 (<x(T)) < 3K and C < £i = (9Ci^) _1 . □ 

Lemma 3.6 Assume that (p, u) satisfying \3.4\ ) with K > as in Lemma 13.51 is a smooth 
solution of \1. 1\) - [T75\) on BL 3 x (0, T]. Then, there exists a positive constant e 2 > 0, depending 
on p and M , such that 

$i(T) + $ 2 (T) < Cl'\ (3.37) 

provided Cq < e 2 . 

Proof By (J33D, and (I3T28D . we find 

o-llVtiH^ 3 < Co- (||/on||L2 + ||Vu|| L 2 + \\p - p s \\ 2 L 6 + ||/9 - />s||l 6 ) 2/3 

<c(c 1/6 + c 2/9 + c 1/9 )<cc 1/9 . 

Thus, using fl37T]), ([Ml), (l3~28t (I3~33l) and (^361) . one derives from fl3~T0D and (I3~TTT) that 
$i(T) + $ 2 (T) 

f T „ „<, r (T) „ ,„ 

<CC + C/ o-\\Vuf L3 dt + C (j||Vn|| 3 3 dt 
Vo-m Jo 



a(T) 

T pa(T) 



L 2 T v || v a\\ L 4j uo -r o / u || v u-11^2 II v "H^e 
<t(T) ' " ./ll 



< CC n + C I (\\Vu\\'i 2 +a 6 \\Vu\\U) dt + C I a\\Vu\\ 3 / 2 \\Vu\\ 3 T /2 dt 
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< cc + cc 1/9 [ a (\\pu\\ L 2 + \\Vu\\ L 2 + ||p - Ps \\ 2 L6 + \\ P - Ps \\ L6 ) 

Jo 

( \ 5/1 

< CC Q + CC 2/3 f J (\\pu\\ 2 L2 + \\Vu\\ 2 L2 + \\p - psWie + \\p - Ps \\l 6 ) dt \ 



< CiC 2/3 < c 1/2 , 

provided Cq < E2 = min{eo 5 £\, C-f 6 }. The proof of Lemma I3~6l is thus complete. □ 
To complete the proof of Proposition 13. 1\ it remains to prove the uniform upper bound of 
the density. To do this, we first need the following refined estimates. 

Lemma 3.7 Assume that (p, u) satisfying \°3.4\ ) with K > as in Lemma \3. 51 is a smooth 
solution of (EZOJ-fOP on M 3 x (0,T]. Then 



sup \\Vu\\ 2 L 2 + / \\p 1/2 u\\ 2 L2 dt < C, (3.38) 

0<t<T Jo 

sup (a\\p 1/2 u\\ 2 L2 ) + [ a\\Vu\\ 2 L2 dt < C, (3.39) 

0<t<T ^ 'JO 

provided Co < £2- 

Proof. We only have to prove (|3.39p since the estimate ()3.38|) is an immediate result of (|3.36p 
and (|3.37p . Choosing m = 1 in (|3,26p and integrating it over (0,T) yields that 

sup ( cr|| / o 1 / 2 -u||^ 2 ) + / crHVnll^dt 
0<t<T ^ ' JO 

<c(c + ^ \\p 1/2 u\\ 2 L2 dt + a\\\7u\\^4dt I 

<C + C a\\Vu\\ L 2\\Vu\\l a dt 
Jo 

<C + C a(\\pu\\ 3 L 2 + \\Vu\\ 3 L2 + \\p- Ps f L 6 + \\p- p s \\ 3 L6 )dt 
Jo 

< C + C sup (V /2 ||p 1/2 tt|| L2 ) , 

where we have used (|3.9p . (|3.33[) . (|3.28p and (|3.38|) . Combining this with Young inequality 
immediately leads to (|3.39p . □ 
We are now ready to prove the uniform upper bound of the density, which is in fact the key 
to extend the local smooth solution to be global and will be proved by modifying the arguments 
in |16 t ll7 1 [2"T] basing on the Zlotnik inequality (cf. Lemma I2.2p and the standard "effective 
viscous flux" F (see (|1.28p ) in a similar form as the one in |lU y i7 t l2"3]. 

Lemma 3.8 Assume that (p, u) satisfying {3.$ with K > as in Lemma 13.51 is a smooth 
solution of i lj.jj) - [775)) on M 3 x (0,T]. Then there exists a positive constant i > 0, depending on 
p and M, such that for all (x, t) € R 3 x [0, T], 

p{x,t) < -p 

provided Cq < e. 
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Proof. Let F be the standard "effective viscous flux" denned by (|1.28p . Then (jl.ip can be 
written in the form: 

D t p = g(p)+b'(t), 

where 

In order to apply Lemma \2.2\ we have to deal with b(t) since lim g(p) = — oo. To do this, 

p— s-+oo 

we first observe from (11.21) and (11.61) that 



AF = div(pu-(p-p s )Vf). (3.40) 
Applying the standard L p -estimate to the elliptic problem ()3.40p gives 

\\VF\\ L 4 < C\\ P U\\ L A + C\\(p - p 8 )Vf\\ L A 

< C\\pu\\]/ 2 %u\\% 4 + C\\p - p s \\ L i2 \\Vf\\ L e (3.41) 

1 3/4 

II 2 



<Co-- % / 8 \\Vu\\% 4 + CC 1 /12 , 



where we have also used (JEBJ), fl33J) and (137391) . It thus follows from (fl~28l) . (12T2D . fl3J|, (13381) 
and (j374T]) that 

\\F\\ L oo < CIIFII^/IIVFH^ 7 

< C (\\Vu\\^ + \\p- p^j) (a-^\\Vu\\% 4 + Co 1 / 12 ) 67 ' 

< Co-V 28 \\Vu\\% 14 + CC] ,l \ (3.42) 
For < tx < t 2 < o-(T) < 1, we deduce from ([33Tj) . (|3T39|) and ([37^2]) that if C < e 2 , then 

/.<t(T) 

< C / ||F|| L oo(it 
J o 

|2 ^3nv7„M|2 



< CC 1/M + C f a" 4 / 7 (a||Vn||i 2 ) 1/4 (a 3 ||Vu||i 2 ; 

JO 



/ r( T ) \ 1/4 / r( T ) \ 1/14 

< CC 1/14 + C ( y a||V«||| 2 dij ij a 3 \\Vu\\l 2 dt ] 



<CC 1/14 + Cd» 4/14 (T)<CC 1/28 . 
Therefore, for any t £ [0,u(T)], one can choose No,Ni in (|2.3p and in (|2.5p as follows: 

n = cc 1/28 , iVi = o, r = p. 

Then, due to the facts that p < p s < p (see Proposition II. ip and 

9(0 -"IT^ 7 "^-"^ ' v ^- r = ^ 
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it follows from (|2.4p that (keeping in mind that < po < P an d /5 > p + 1) 

sup Hp^IUoc < max{p, p} + N < p + CC Q 1/28 < -p, (3.43) 

Q<t<a(T) 1 

provided the initial energy Co is chosen to be such that 

\ 28 
2CJ ' 

For any a(T) <h<t 2 <T, we have from (|3Ji, (l337j) and (|33Hj) that 



C < min{e 2 ,e3} with e 3 = ( 



\h{t2)-h(t X )\ 

< C [ 2 \\F\\ L oodt < C [ 2 ||F||^ 2 7 ||VF||^ 4 7 cft 



< cc 1/14 (t 2 - tl ) + c r \\ P ^% u \\u\\% u dt 

Ju 



< ( CCl /li + (t 2 -h) + C £ a 3 (\\p 1/2 u\\h + Wp4» ) d t 



2{2p + A) 

A 



where in the last inequality we have chosen Co to be such that 
Cq < min{e2, £3, £4} with £4 -- 



A - u 



2C{2p + \) / 

Therefore, for any t 6 [a(T),T], we can choose Nq and N\ in (12.3P and in (12.5|) as follows: 

iV = CC 1/2 , Ar 1 = _^, e* = P + l. 
Then it is easy to check that for any ^ > £*, 

<?(6<-77^t(£ 7 -^)<-^i 1 



2// + A r ' - 2/i + A 

and hence, one has from (|2.4p and (|3.43[) that 

sup <max\h,p + l\ +N < ls + CC 1/2 < U, (3.44) 

a{T)<t<T I 2 J 2 4 

provided the initial energy Co satisfies 

/ ~ \ 2 

C < e = m.m{e2,s 3 ,S4,e 5 } with e 5 = ( — 
The combination of (|3.43p with (|3.44j) proves Lemma 13.81 □ 
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4 Time-dependent higher-order estimates 

In this section, we prove the higher-order estimates of the smooth solution (p, u) to (jl.lD - 
(jl.5p . which are needed for the existence of classical solutions. Prom now on, we always assume 
that the initial energy Co satisfies (|3.6p and that the conditions of Theorem 11.11 are satisfied. 
We also denote by C the various positive constants which may depend on po,uq, ||<?||z,2,/, A, p, 
A, 7, p, p, p, M, and T as well. 

We begin with the L 2 -estimate on the material derivative of the velocity. 

Lemma 4.1 For any given T > 0, there exists a positive constant C(T) such that 

sup \\p 1/2 u\\ 2 L 2 + [ \\Vu\\ 2 L2 dt<C(T). (4.1) 

0<t<T ' Jo 

Proof. Choosing m = in (|3.26[) and integrating it over (0, T), we have 



sup / p\u\ 2 dx — / p\ii\ 2 (x, 0)dx + / ||V-u||^2<ii 

<t<TJ J JO 

<cJ o (l|V«||ia + \\p 1/2 u\\ 2 L2 ) dt + C J o H V 

<C + C [ \\Vu\\ L 2\\Vuf L edt 
Jo 

r T 

<C + C \\puf L2 dt 
Jo 



u\\ 4 L 4dt 



<C + C sup \\p 1/2 u\\ L 2, 

0<t<T 

where we have used (|3.7p , (|3.9p , (|3.28p and (|3.38p . Combining this with Young inequality gives 
(|4.ip . since the compatibility condition (|1.14p implies that ^/pu\t=o = ^/po{^ f — g) € L 2 is well 
defined. □ 
Next, similar to that in [17, 18J, we utilize the Beale-Kato-Majda-type inequality (see Lemma 
2.3b to prove the important estimates on the gradients of (p,u). 



Lemma 4.2 There exists a positive constant C = C(T) such that 

[\\Vp\\vmfi + \\Vu\\w)+ I 

0<t<T 



sup {.\\Vp\\ L 2 nL6 + [|V«|| H i)+ / \\Vu\\ L ™dt < C(T). 

)<t<T Jo 



(4.2) 



Proof. Since C = —pA — (p + X)S/div is a strong elliptic operator (see [3] for instance), applying 
the L p -estimate of elliptic system to (|1.2|) gives that for 2 < p < 6, 

||V 2 n|| LP < C (\\pu\\ LP + \\Vp\\ LP + \\Vf\\ LP ) . (4.3) 

By integration by parts, we easily derive from (jl.ip that 

^||Vp|| LP < C\\Vu\\ LOO \\X7p\\ LP + C||V 2 n|| iP 

< C (1 + ||V«||l») ||Vp||iP + C (1 + ||^u||lp) , (4.4) 
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where (14. 3p was used in the second inequality. Using (|3.38p . (|4.3|) and (12. ip . we deduce from 
(1231) that 



|| Vm||l°c < C (||divu||ioo + ||curln||Loo) In (e + ||p?i||L6 + ((VpH^) + C 
<C + C (Hdivullioo + ||curln|| L oo) In (e + ||Vn|| L 2) 

+C(||divu|| L oo + ||curln|| L oo)ln(e + ||Vp|| LB ). (4.5) 

Set 

$(t) = e + ||Vp|| L6 , V(t) = 1 + ([|div«||r« + ||curl«||ioo) In (e + ||Vtt|| L2 ) + ||Vtt|| L 2. 
Then it follows from (|4.4p with p = 6 and (|4.5p that 

< C*(i)$(i) In 
which, together with the fact that <&(t) > 1, implies 

— ln$(t) < C*(t)ln$(t). (4.6) 

Applying the standard L p -estimate to (13.311) and (|3.32p yields that 

\\VF\\ L e + ||V( / 97 1 curlu)|| i 6 < C {\\fm\\ L e + ||Gi|| L e + ||G 2 || L 6 + || VuVpJ 1 \\ L e) 

<C(\\Vu\\ L 2 + ||V/£> a |Uoo||V«|Ue + ||VpJ £6 ) 
<C(||Vn|| L2 + l), 

where we have also used CL11, §SM, (1338]) . (I3T38D . (I4TT]) and Lemma Thus, 

I ^{t)dt<C + C f {\\Vu\\ 2 L 2 + ||div«|||oo + HcurluHloc) dt 
Jo Jo 

<C + C [ (||P(^)-P(^)||i- + ||^||icc + ||^ 1 cuxl«||ioc)£tt 
Jo 

<C + C [ (||F||2 2 + ||curlw||2 2 + ||VF||^ 6 + HV^WmUs) 
Jo 

<C + C f \\Vu\\ 2 L2 dt < C, (4.7) 
Jo 

and consequently, it follows from (I4.6P and Gronwall's inequality that 

sup \\Vp(t)\\ L6 < sup $(t) < C. (4.8) 

0<t<T 0<t<T 



As a result of (|43|) . (|1T7|) and dMD, we obtain 

/ ||Vn|| L oodt < C + C [ ^{t)dt<C. 
Jo Jo 

Using (14.ip and (14.9p . we infer from (14. 4p with p = 2 and Gronwall's inequality that 



(4.9) 



sup ||Vp(i)|| L2 <C, 

0<t<T 

which, together with (|4.ip and (|4.3p . also implies that ||Vn||^i < C. The proof of Lemma 14.21 is 
therefore completed. □ 
By virtue of Lemmas 14.11 and 14. 2\ we easily obtain 
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Lemma 4.3 For any given T > 0, it holds that 

cT 



sup \\p 1/2 u t \\ 2 L2 + f \\Vu t \\ 2 L2 dt < C(T), (4.10) 

0<t<T JO 

sup (\\Vp\\ m + \\VP\\ H i) + f T \\V 2 u\\ 2 m dt < C(T), (4.11) 

0<t<T 70 

sup (\\p t \\ H i + \\P t \\ H i) + f (||p tt || 2 2 + \\P tt \\ 2 L 2) dt < C(T). (4.12) 

0<t<T JO 

Proof. First, (|4.1U|) follows directly from Lemmas 14.11 and 14.21 
Next, we prove (|4,lip . Since P{p) = Ap 1 satisfies 

P t + u- VP + jPdivu = 0, (4.13) 

from which and (jl.ip we have by f|4. 2j) and direct computations that 

|(l|V 2 ,||i 2 + ||V 2 P||i 2 ) 

< C||V«[|r« (l|V 2 p||i 2 + ||V 2 P||i 2 ) +C\\V 3 u\\ L 2 {\\V 2 p\\ L 2 + ||V 2 P|| i2 ) 
+C\\V 2 u\\ L3 (\\V P \\ L6 + ||VP|| L «) (||V 2 p|U 2 + ||V 2 P|| i2 ) 

< C(l + ||V«||l<») (IIVVH^ + ||VP|| 2 2 ) + C||V\||^i. (4.14) 

Using (|4.2|) . (|4.1U|) and Lemma [2ZU we deduce from (jl.2j) and the standard L 2 -estimate of elliptic 
system that 

\\V 2 u\\ m < C {\\ P u t \\ m + \\pu ■ Vu\\ m + \\VP\\m + ||pV/|| H i) 
< C (1 + ||Vn t || L2 + HVplUsKlUe + ||V 2 P|| i2 ) 

<C(l + ||Vn t || L2 + ||V 2 P|| L2 ). (4.15) 

Thus, putting (|4.15p into (|4.14p and using Gronwall inequality, we immediately arrive at ()4. 
since it holds that [|Vit||i,«. + ||Vut|| 2 2 G F X ($,T) due to §£2$) and (l4TT0]l . 

Finally, we prove (|4TT2|h Thanks to (g2]) and ([ITTT]) . it follows from ([LTD and (|iT3|) that 

\\p t \\ L 2 + ||P|| i2 < C||u|| £ « (\\Vp\\ L 2 + \\VP\\ L 2) + C\\Vu\\ L 2 < C, (4.16) 

HVptlUa + ||VP t || i2 < C||V 2 n|| i2 + C||u||i=o {\\V 2 p\\ L 2 + ||V 2 P|| i2 ) 

+C||Vn|| L3 (\\Vp\\ L6 + \\VP\\ L e) < C, (4.17) 

where Sobolev inequalities were used to get that ||ii||z,°° + ||V«||x,3 < C||Vii||#i < C. Moreover, 
since (|4.13p implies 

Ptt + Ut ■ VP + u • VP + 7P t divn + 7 Pdivn t = 0, 
one obtains after using (|4T2]> . (liTTOt (I4TTTT) . (|4TT6"]> . and (|4TT71) that 

/ ||Pt|li 2 ^<C / (|K||L6||VP|U3 + ||VPt|U2 + ||Vu||i=x>||P t || i2 + ||V« t || L 2) 2 ^ 

Jo Jo 

<C + C [ (\\Vu\\ 2 H2 + \\Vu t \\ 2 L2 )dt<C (4.18) 
Jo 
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In the same way, one also has Hp^H^a G L 2 (0, T). So, combining this with (|4.16|) - (|4,18|) completes 
the proof of (0112$. □ 
In order to prove the solution obtained is indeed a classical one on the time-interval [r, T] 
for any 0<r<T<oo, we need some further estimates on the higher order derivatives of (p, u). 
However, due to the weaker compatibility condition (jl,14p (cf. (|1.20p . (jl.2ip ). the methods used 
in [T7] cannot be applied any more. To overcome this difficulty, we need the following initial-layer 
analysis. 



Lemma 4.4 Let a = min{l,t}. Then it holds for any given T > that 

rT 

o<t<T "" "" "~ ' Jo 

Proof. First, differentiating (II. 2p with respect to t gives 



sup a {\\Vuf H2 + \\Vu t \\y) + [ T * (\\p l/2 u tt \\ 2 L2 + \\Vu t \\ 2 Hl ) dt < C(T). (4.19) 

)<t<T Jo 7 



uAu t + (p + X)Vdivu t 

= pu u - PtVf + ptu t + ptu ■ Vu + put -Vu + pu- Vu t + VPt- 



(4.20) 



ld_ 

Ydt 



Multiplying (|4.2U|) by Utt and integrating the resulting equality over M 3 yield 

[p|Vii t | 2 + (p + A)(divti() 2 ] dx + J p\u u \ 2 dx 

j to V/ - - „,„ Vu — put Vu — pu Vu, - VP,) ■ uudx 
d /' 1 

-~Pt\ut\ 2 + (ptVf - p t u -Vu -VP t ) ■ u t dx 
-J(pttVf-Puu.Vu-ptu t .Vu- P tu.Vut).utdx 
+ \ j Pu\u t \ 2 dx - j Ptt^utdx - j {pUt -Vu + pu- Vu t ) ■ u u dx 



Then, we estimate each term on the right-hand side of (|4.2ip . Using (jl.ip and integrating 
by parts, we see from P~2j) and (^T0]) - (|4TT2]) that 

lo = / [-PU ■ Vut -ut + ptiVf-u- Vu) • u t + Pt^ut] dx 

<C\\u\\ L o B \\p 1 / 2 U t \\ L 4Vu t \\ L 2 +C||p t || i2 ||V/|| i 3|K|| L 6 

+C||/ J i||L 2 ll u llL°°ll V ' u llL 3 ll u illL 6 + C||P|lL 2 ll Vu tllL 2 

<^||Vn t ||| 2 + C. 
Using (|4T2|) and (|4.10p - (|4.12p again, one infers that 

h < cWpttWtfWv fWtf IKIUs + c||pttllL2|MlL° o l|Vu|| L 3||'ut|| i 6 

+C||PtlU 2 ll«t|ll,6||V«|| i 6 + CllptlliellltllisllVlttll^lllitll^e 
< C\\ptt\\ L A\Vut\\ L 2 + C\\Vu t \\ 2 L2 < C {\\ P tt\\ 2 L2 + \\Vut\\ 2 L2 ) . 
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By (jl.ip and integration by parts, we have from (|4.10p and (14.12H that 



h = J (P u )t ■ Vut ■ u t dx = J {p t u ■ Vu t ■ u t + pu t ■ Vut -Utjdx 

< C , || / 9t|| L 3||n|| L oo||Vn t || L2 ||'u 4 || L 6 + C||pii i || i 3||Vnt|| Z 2||u t || Z 6 

< C\\Vut\\h + CIIPWtll^Kll^llVutH^ < C (1 + HVutHiO . 
Obviously, I 3 < C (\\P t t\\ 2 L2 + ||V«t||| 2 )- Finally, it follows from (j4.2|) and Lemma [2TT1 that 

h < C (\\u t \\ L (i\\Vu\\ L 3 + [|u[|i<x»[|Vut||£,2) \\p 1/2 uu\\l2 
<l\\p 1/2 uu\\h+C\\Vu t \\ 2 L2 . 

Thus, putting the estimates of It into (|4.2ip . multiplying the resulting inequality by a(t) 
and integrating it over [0, T], we infer from (j4.10|) . ()4.12p and Gronwall's inequality that 



0<t<T 

which, together with (|4.1ip and (|4.15p . gives 



sup (<7||V«t|||a) + f a\\p l l 2 u tt \\ 2 L2 dt < C, 

<t<T JO 



sup a (\\Vuf H2 + \\Vu t \\ 2 L2 ) + [ a\\p^ 2 u tt \\ 2 L2 dt < C. (4.22) 

0<t<T JO 



Finally, applying the standard L 2 -estimate to the elliptic system (|4.20p together with (|4.2p 
and (|4. 12|) gives 

||V 2 n t || i2 < C\\pAu t + (p + A)Vdivu t || L 2 

< C\\pu tt - ptVf + Ptu t + ptu • Vn + pu t • Vu + pu- Vu t + VP t \\ L 2 

< C (\\p ll2 U t t\\ L i + \\Pt\\l? + HPtllLslKHie + ll/^IUsNliHI^HIz, 6 
+ IKIIl,8||Vu|| i 3 + ||^||L°°||Vn t || L 2 + ||VP t || L2 ) 

<c(l + 1^11^ + 1^/2^11^, (4.23) 

from which and (|4.22p it follows that 

fT 

/ a\\V 2 u t \\ 2 L 2dt < C. 
Jo 

This, together with (j4T22j) . finishes the proof of (|4TT9|h □ 
The next lemma is concerned with the VF 1,9 -estimate (q £ (3, 6)) on the gradients of density 
and pressure, which in particular indicates the Holder continuity of (Vp, VP). 

Lemma 4.5 Let q E (3,6) be as in Theorem \l.l[ For any given T > 0, it holds that 

sup (\\Vp\\ m , g + ||VP|| W1 „) + [ T (|M|$i„ + ||V 2 u||^ M ) dt < C(T), (4.24) 

0<t<T Jo 

where 

po 4 (9q - 6)/(10 9 - 12) G (1, 4q/(5q - 6)). (4.25) 
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Proof. Applying the differential operator V 2 to both sides of (|4.13p . multiplying the resulting 
equations by q\ V 2 P(p)\ q ~ 2 V 2 P(p), and integrating it by parts over M 3 , one deduces from (|4.2p . 
(jITTj) and Lemma O that 

j t \\^ 2p \\h < c {\\vu\\ L ~\\v 2 p\\ Lq + ||vp||l-||v 2 «||l, + \\v 2 u\\w^) llv'Piil; 1 
< c (i + \\Vu\\ H .) (i + ||v 2 p||1,) + ciiv^ii^nv^iil; 1 . 

The similar estimate also holds for ||V 2 /)||l9. Therefore, 

f t (l|V 2 p||| 9 + || V 2 P||1 9 ) < C (1 + || V«||flO (1 + ||V 2 p||l 9 + || V 2 P||l g ) 

+C||V 2 u|| w x, 9 (llVVHl; 1 + ||V 2 P||l; 1 ) . (4.26) 



Applying the standard VK 1,p -estimate to the elliptic system (jl.2p yields that 

||V 2 U|| W 1, 9 < C (\\ P U t \\ w i, q + \\fM ■ Vu\\ w i, q + || VPHwrl,, + \\pVf\\ w i, q ) 

< C(l + \\ut\\ W i, q + ||Vp|| M |KIU- + UVullwri,, 

+ ||Vp|Uoo||Vu|U a + ||Vu|Uoo||Vu|U« + HVPHhti,,) 

< C (1 + ||V«||fla + \\u t \\ w i, q + \\Vp\\ m , q + \\VP\\ w i, q ) , (4.27) 

where we have also used (|4.2p . (|4.1ip and Lemma l2.1i Putting (|4.27p into (|4.26p gives 

j t (llVVlli, + l|V 2 Pf L9 ) < C (1 + HVuHfl, + Hutllnri.,) (1 + \\V 2 P \\ q Lq + \\V 2 P\\l q ) . (4.28) 
Using d2H, (J221), (USD and (|PTi|l . we find that 

r htwhdt < c r \\u t f H1 dt <c r (m*, + iiv^m,) * 

JO JO JO 

(||p 1/2 n t |li 2 + ||V« t || 2 2 ) eft < C. (4.29) 



Note that 4q/(5q — 6) G (1,4/3) for g E (3,6). So, for po as in (|4.25p . we obtain by Lemma [2TTT 
Holder inequality and (|4.19p that 

\\VutWl\dt < C f T UVutll^-^/^UVtitll^-^^cft 

JO 

< c [ T a-^ 2 {a\\Vu t \\ 2 L2 ) M(i - qmq) (a\\Vu t \\ 2 Hl ) po{3q ~ m4q) dt 

Po(6-g)/(4g) 







2 xp (3g-6)/(4 9 ) 

T \ (4g-p (3g-6))/(4g) / T \ Po(3g-6)/(4g) 



< C sup a||Vu*|| 2 2 / a^' 2 HV^) 90 ^'™ dt 

\0<t<T I JO 



(pi \V*l—P0\,'Xi—vjJ/V*d) / pT \i 

I a -2 PQg /(4 g - P0 (3 g -6)) dt \ I a\\Vu t \\ 2 Hl dt) 

< C, (4.30) 
since < 2p q/(4q - p Q (3q - 6)) < 1 and < p (3q - 6)/(4g) < 1. 
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The combination of with (OD1) shows that for q £ (3, 6), 



IKH^dt^C. (4.31) 
Thus, by ([4, lip . (|4,31|) and Gronwall's inequality, one sees from (|4.28|) that 

SUp (IIVpllw-l,, + HVPHiyl,,) < C, 
0<t<T 

which, combining with (|4.1ip . (j4.27|) . and ()4.31 j) . finishes the proof of Lemma 14.51 □ 
Finally, we still need the following lemma, which implies that Ut and V 2 u are Holder 
continuous away from t = 0. 

Lemma 4.6 For any given T > 0, it holds that 

sup a (\\p 1/2 u tt \\ L 2 + \\V 2 u t \\ L 2 + \\V 2 u\\ w ,, q ) + / a 2 \\Vu tt \\ 2 L2 dt < C(T). (4.32) 

0<t<T v ' Jo 

Proof. Differentiating (|4.20p with respect to t gives 

puttt + pu ■ Vutt - pAuu - (p + A)Vdivu 4t 
= 2div( / ou)ii tt + div( / on) t u i - 2(pu) t ■ Vu t 

-(pttu + 2p t u t ) ■ Vu - puu ■ V-u - VP tt + PttV/, 

which, multiplied by u tt in L 2 and integrated by parts over R 3 , yields 
~ J p\u tt \ 2 dx + J (p\Vu tt \ 2 + (p + A)|divu tt | 2 ) dx 



pu ■ Vu tt ■ u tt dx - J (pu) t • (V(ut • u tt ) + 2Vu t ■ u tt ) dx 
- j (p tt u + 2p t ut) ■ Vu • uttdx - J putt • Vn • u tt dx 

+ fPtt^u tt dx + JpttVf- u a dx ±± J, (4.33) 

i=i 

The right-hand side of (|4.33p can be estimated term by term as follows, using Lemma |2. H 
Cauchy-Schwarz inequality and the estimates obtained. 

Ji < \\p 1/2 utt\\ L 2\\Vu t t\\ L 2 < S\\X7u t t\\ 2 L 2 +C{5)\\p 1 l 2 utt\\ 2 L 2, 
h < C(\\pu t \\ L 3 + WptWtf) (l|Vn t || L2 ||«tt|| L 6 + |K||i,6||Vutt|| L 2) 
< C(l + ||pn t ||^ 2 ||Vn t ||^ 2 ) \\Vuth4Vutth* 



<6\\Vuttf L 2 + C(S)(l + \\Vut\\l 2 ), 



'h<C {\\Ptt\\LA\ U \\L^\\^ U\\ L 3 + ||/0t||i3||'Ut|| i 6||V«|| i 3) ||«tt||i6 

<J||Vn tt || 2 2 +C(5)(||^|| 2 2 + ||Vn t || 2 2 ), 
Ja < C||p 1/2 n tt || L2 ||V U || L 3||n ti || L 6 < 5\\Vutt\\ 2 L 2+C(5)\\p l ' 2 u t t\\ 2 L 2 
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and 

J 5 + h <C(\\P U \\ L 2 + \\p u \\ L i) (\\Vu tt \\ L 2 + \\Vf\\ L3 \\u tt \\ L6 ) 
<5\\Vuu\\h + C(8)(\\Pu\\h + \\pu\\h). 

Putting the estimates of J\ , . . . , Jq into (|4.33j) , multiplying it by a 2 and integrating the 
resulting relation over (0,T), we have by choosing 5 > small enough that 



sup (<T 2 \\p 1/2 u tt \\ 2 L2 ) + [ a 2 \\Vu tt \\ 2 L2 dt 

0<t<T v 7 JO 

<c + c£ (<j\\ P ll2 u tt \\h + ° 2 (\\Pu\\h + \\p*\\h + WVutWhj) 

<C + C sup (a^WVutWlA f T \\Vut\\ 2 L 2dt < C, (4.34) 

0<t<T v ' Jo 

where we have used (ETlOD . (j4TT2l) and ()4TT9l) . 

As a result of P~34"j) . (|4TT9|) and ([^723]) . we also see that 

^||v 2 n t || L2 ^^(i + iiv^ii^ + Hp 1 / 2 ^!!^) <c, vte[o,r]. (4.35) 

and thus, it follows from (ET2T1) . (1^191) and ()4T24"D that for any t £ [0,T], 

0-||V 2 -u||jyi, 9 < Ca(l + \\Vu\\ H 2 + ll^ll^/i.g + ||Vp||^i,9 + ||VP|| w i, ? ) 

<C(l + a\\u t \\ H 2) <C, (4.36) 

where we have used (f277|) . (f379|) . (|4.10p and (|4.19p to get that <7||ii t || L 2 < C. Lemma HUl now 
readily follows from (Q4]l - (l436]l . □ 

5 Proofs of Theorems 11.11 and 11.21 

With all the a priori estimates at hand, we are now ready to prove our main results. First of 
all, we prove the global well-posedness of classical solutions to the problem (|1.1|) - ([1.5|) provided 
the initial density is strictly away from vacuum and the initial energy is small. 

Proposition 5.1 Let Co be the initial energy defined by U.10\) . For given numbers M > (not 
necessarily small) and p > 0, assume that po,uo and f satisfy 

'feH 3 , (p -Poo,u ) G H 3 , 

< inf po(x) < sup po(x) < p, \\Vuo\\l 2 < M, (5-1) 
. Co < e, 

where e > is the same one as in Proposition \3.1[ Then for any < T < oo, there exists 
a unique classical solution (p,u) of ( II. 1\) - 1775)) on M 3 x (0,T] satisfying 112. 9\) . 112. 1 0\) with Tq 
replaced by any T > 0. Moreover, all the uniform-in-time estimates in Sect. 3 hold for (p,u). 



dt 
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Proof. The standard local existence result (i.e. Lemma 12. 5p shows that the Cauchy problem 
(|l.ip - (|1.5p has a unique local classical solution (p,u), defined up to a positive To which may 
depend on inf po(x) and satisfying (12.9p . (I2.10p . 

In view of (|5.ip and the definitions of 3>i(T) (i = 1, 2, 3), we know that 

$!(()) = $ 2 (0) = 0, $3(0) < M and p < p. 

Thus, there exists a positive T\ G (0, To] such that (13. 4p holds for T = T\. 
Set 

T* = sup{T I jSap holds} . (5.2) 

It is clear that T* > T\ > 0. 
We claim that 

T* = 00. (5.3) 

Otherwise, T* < 00. Then it follows from Proposition 13.11 that (|3.5p holds for any < T < T*, 
and furthermore, the estimates in Lemmas I4.1H4.6I are also valid for all < T < T*. 

In the next, for the sake of simplicity, we denote by C the positive constant which may 
depend on the lower bound of the initial density (i.e. inf po(x)) and T*. We shall prove that 
there exists a positive constant C > such that for any T G (0,T*), 

sup ||(p-/3oo,w)||jf3 < (7. (5.4) 

0<t<T 

This in particular implies that 

|| (p - poo, u)(-, T*)\\h3 < 00, inf p(x,T*)>0. 

So, Lemma [23j together with (|3.5p . yields that there exists some TU > T* such that (|3.4p holds 
for T = T**, which contradicts (I5.2p . Hence, (I5.3P holds. Thus, (p,«) is in fact the unique 
classical solution of pH]) - (fT5]) on E 3 x (0, T] for all < T < T* = 00. 

We are now in a position of proving (I5.4p . To do so, we first note that under the conditions 
of f)5. 1 [) 1 and ()5.ip 9. it holds that (keeping in mind that po > 0) 

ut(; 0) = V/ - n • Vu + p 1 (pAn + (p + A)Vdivn - P(p )) G H 1 . (5.5) 

Thus, similar to the proof of f|4.22j) . by (|5.5p one easily deduces that 

rT 

0<t<T ' JO 

Similar to the derivation of (14.231) . one gets by using (j5.6[) that 



sup (\\Vuf H2 + ||Vu t ||| 2 ) + f \\p^ 2 u tt \\l 2 dt < C. (5.6) 

Kt<T JO 



£ \\V 2 u t \\ 2 L2 dt < c£ (l + ||Vn t ||2 2 + iipVV,^^ dt 



r 



<C + C \\p 1/2 u tt \\ 2 L2 dt<C. (5.7) 







Moreover, using the standard H 2 -theory of elliptic system, one deduces from Lemmas 12.11 14.11 - 
^3] and ([511 tnat 

||V 2 n|| H 2 < C\\pAu + (p + A)Vdivn|| jff 2 



29 



< C\\pu t + pu-Vu + VP - pVf\\ H 2 

<C(\\pu t \\ L 2 + IIVpIIhiIKHloc + \\Vp\\ LS \\Vu t \\ L e + \\Vu t \\ H i 

+ \\PU\\ H 2\\VU\\ H 2 +\\VP\\ H 2 + \\Vf\\ H 2 + ||Vp|| H i||V/|| H2 ) 

<C(1+ \\V 3 p\\ L 2 + \\V 3 P\\ L2 + \\Vu t \\ m ) . 



(5.8) 



Applying the differential operator V 3 to both sides of (jl.ip and (I4.13p . and multiplying the 
resulting equations by V 3 p and V 3 P(p) in L 2 , respectively, one easily obtains after integrating 
by parts and using Lemmas 12.11 14.1H4.3[ (|5.6p that 



d_ 

di 



^ p\\h + ll v P\\h) <c{i + \\v 2 uf H 2 + \\v 3 P \\ 2 L 2 + ||v 3 p|| 2 2 ) . 



Combining (|5.7p . (|5.8p and (|5.9p . we conclude from Gronwall's inequality that 

sup (||v 3 p||| 2 + ||v 3 P||£ 2 ) < c. 

0<t<T 



(5.9) 



(5.10) 



So, Lemmas 14. 1H4.3| together with (|5.6p and (|5.10|) . lead to Q5.4|) . The proof of Proposition 15.11 
is therefore complete. □ 

Proof of Theorem \l.l[ Let (pq,uq) be the initial data satisfying the conditions (j!.12|) — (j!.14p 
in Theorem II .11 Assume that the initial energy Cq satisfies 



Co < e = e/2, 

where e > is the same one as in Proposition 13.11 

To apply Proposition 15. 1\ we define the smooth approximate data as follows: 



(5.11) 



Po 



6,r, _ 38 * Po + V 



5,7] 



35*9, f ,v = 35 * f, 



(5.12) 



1 + 7? ' 

where < S, r] < 1, jg(x) is the standard mollifier with width 5, and " * " denotes the usual 
convolution operator. As that in [2], let « ' ?? be the unique solution to the elliptic problem: 



-pAu 6 " -(p + A)VdivuJ" = -VP(p^) + [p S ' v 
^ Uq ,ti (x) —t-0 as \x\ — > oo. 



1/2 



(5.13) 



Also let p S s ,v = p S s ' r] (x) be the unique solution of (|1.6|) with / being replaced by f 6 ' v . Obviously 
pi' 71 satisfies (|1.8p . (|1.9p . Moreover, it is easy to check that 

' < < pj" < p < oo, 

(p^-Pocg^^J^^pi^eH^, 

lim ( \\p S Q ,v - p \\ H 2 nW 2, q + \\u ,v - u \\ H 2 ) =0, 

5,r)— >0 \ / 

lim (\\g s * - g\\ L 2 + IK/ 5 '" - f,p*« - Ps )\\ H 2 nW 2, q 
The initial norm of the mollified data (pg 7 ', it^) now reads 

cf " = / (go*") + ^'>o'T) dx, 



(5.14) 



0. 
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where G(-) is the function defined in (II. lip with p s replaced by pt' 71 . 
By (|5.14|) . we have 



lim lim C, 



5,T) 





C . 



So, there exist positive constants 770 G (0, 1) and £0(77) such that 

< Co + e/2 < i, 



(5.15) 



provided that 

< i) < rjo and < 5 < 5q(tj). (5.16) 

Let 5,7] satisfy f[5.16]> . Then it follows from (|5.15p and Proposition 15.11 that there exists a 
unique smooth solution (p s,v , u s,v ) of (jl.ip . (jl.2p with the mollified data p ,ri u n ,ri , f 6 ' 71 and g 5,7 ' 
on M 3 x (0, T] for all T > 0. Moreover, Proposition 13.11 and Lemmas I4.1H4.61 independent of 8 
and 77, hold for {p 5,n , 11 s ' 71 ). 

Now, passing to the limit first 5 — > 0, then 77 — > 0, we deduce from standard arguments that 
there exists a unique solution (p, u) to the origin problem fjl . lj) — (jl .5j) on M 3 x (0, T] for all T > 
satisfying 

' < p{x, t) < 2p for all x G M 3 , t > 0, 
( /0 -p oo ,P(p)-P( /0oo ))GL° o (0,T;F 2 n^), 

(5.17) 

?j G L°°(0, T; H 2 ) n L 2 (0, T; # 3 ) n L°°(r, T; F 3 n VF 3 ><?), 

k «t G L^i^T-H^nH^^T-H 1 ), 

for any < r < T < 00. Note that, the uniqueness of (/?, u) satisfying ()5. 17[) can be proved in a 
standard way as that in [2]. 

To complete the proof of the first part of Theorem 11.11 we still need to show that 



(P 



,P(p)-P{ Poo )) eC([0,T};H 2 nW 2 ^), uGC([0,T];H 2 ). 



(5.18) 



The proof of (|5.18|) is similar to that in [13], and we sketch it here for completeness. First, by 
virtue of fO|l . (|4TTU1) - ([4TT2D and (g^U), it is easy to get that 



f (p - Poo ,P(p) - P(poo)) G C([0, T]; P 1 n PF 1 ^), 
(p-poo,P(p)-P(Poo)) G C([0,T];tf 2 nW^-weak), 
ueC([0,T];H 1 nW 1 'i). 

Denote Aj - d% with i,j = 1,2,3. Note that it holds in T>'(R 3 x (0,T)) that 

dtDijp + div(iiZ)yp) = —div(pDiju) — div(dip • djU + djp • diu). 
Let j^a?) be the standard mollifying kernel with width u, and set p v = p* j v . Then, 
d t Dijp" + dw(uDijp") = -div(pDiju) * j v - div(dip ■ djU + djp ■ d{u) * j v + R u , 
where R v = div(uDijp u ) — div(uDijp) * j v satisfies (cf. [231 Lemma 2.3]) 



(5.19) 



T 



\Rv\\L 2 r\L<idt ^ C 



li/ll po II 11 -n\\ po rii < (1 



(5.20) 



(5.21) 
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with po > 1 as in (|4.25p . 

Multiplying (|5.20p by q\Dijp u \ q ~ 2 Dijp u and integrating by parts over M 3 , we obtain 



d „ ri 



v\\q 



dt 



«jr WLi 



-(q—1) J \Dijp u \ q d\\udx — q j div(pDiju) * j u \Dijp u \ q 2 Dijp u dx 



-q j (div(d i p-d j u + d j p-d i u)*j u )\D ij p u \ q 2 D i!j p v dx 
+q J R v \D ij p v \ q ~ 2 D ij p u dx, 
which, combining with (|4.11j) . (|4.24p and (|5.2ip . yields 



sup iivviil+ r 

0<t<T ' Jo 



,,'2^119 

dt" ' 



VZ i 



Po 

(ft 



<C + C [ T (\\VuQ 2tq + ll^ll^) < C. 

JO 



This, together with Ascoli-Arzela theorem, gives 

\\X7 2 p u (;t)\\ Lq ^\\V 2 p(;t)\\ Lg in C([0,T]), as 1/ -+ 0, 
which particularly implies 

||V 2 p(-,*)IU« GC([0,r]). (5.22) 

Similarly, one can also obtain that 

\\V 2 p(;t)\\ L ,eC([0,T}). (5.23) 
So, it readily follows from (ET22]) . ([£23]) and (EUQlb that 

v 2 pe C([o,r] ; L 2 nL 9 ). (5.24) 

In the exactly same way, we also have 

V 2 P(p) eC([0,T};L 2 DL q ). (5.25) 

It follows from (fL2|) that 

(put)t = (/iAu t + (jjt + A)Vdivn t ) + #V/ - VP* - (pu • V«) t . 
Thus, it is easily seen from P~2]h P~T0j) and P~T2]) that 

||(put)t||ff-x <C(l + ||Vnt|U 2 ), 

and hence, 

(jmt)tZL 2 {Q,T;H- x ). 
This, combining with the fact that pu t £ L 2 (0,T;F 1 ) due to (|4.10p and (|4.1ip . leads to 

W £C([0,T];L 2 ). (5.26) 

Thus, by virtue of (|5.19p . (|5.24p and (|5.25p . one easily gets 

u 6 C([0, T]; H 2 ), (5.27) 
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since it holds that 



liAu + (p + A) Vdivu = pu t + pu-Vu + VP(p) - pVf. 



Combining (157T91) . (Oil) . (f5T25|) . and (IQ71) . we finish the proof of ([5TT8j) . 

To complete the proof of Theorem 11.11 it remains to prove (|1.18p . To do this, we first 
deduce in a manner similar to the derivation of (|3.34|) that 



s IIp-*IL- 



dt 



< 



(5.28) 



/oo 
{\\Vu\\ 2 L2 + \\p- Ps \\t 4 + \\F\\ 4 L4 )dt<C, 

where we have used (|3.7p and (|3.33p . Thus, it follows from (|3.33p and (|5.28p that 

||p — p s \\la — > as t — > oo. 
This, together with (|3.9p . Proposition [I7T] and the interpolation inequality, immediately gives 

\\p — Poo\\lp — > as t — > oo, VpG(2,oo). (5.29) 
To study the large-time behavior of the velocity, we set 

M(t)^||V< 2 + ^||div< 2 . 

Then, multiplying ()3.12p by u in L 2 and integrating by parts, we obtain 

\M'{t)\ <C j p\ii\ 2 dx + C||Vn||| 3 + C||Vn|| L 2. 

Here, we have used (|3.14|) . (|3.15|) and the following estimates: 



J u-V{P(p)-P( Ps ))dx 



<C\\p- Ps \\ L 4Vu\\ L 2 <C\\Vu\\ L 2 



and 



J(fi~ Ps) 



u ■ Vfdx 



u\\ L 2 ) dt 



< C\\p-p s \\ L 2\\u\\ L6 \\Vf\\ L 3 < C\\Vu\\ L 2 
Thus, it follows from (I3TH31 . (fBTBTD . (|3T39l and ([5301) that 

/oo roo . 

\M'{t)\ 2 dt <C J [\\p 1/2 u\\l 2 + ||Vu||^ 2 ||Vu||^4 + || V 

< C J (||p 1/2 u||£ 2 + ||Vu||^4 + ||Vw||| 2 ) dt < c. 
Due to (|3.38p and (|3.7p . one has 

/■oo poo 

/ \M(t)\ 2 dt <C sup \\S7u\\ 2 L2 \\Vu\\ 2 L2 dt < C, 
Ji t>i Ji 

from which and (|5.3ip . we know that 

||Vw(t)|| L 2 — > as t — > oo. 



(5.30) 



(5.31) 



(5.32) 
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As a result, we also have 

j p 1 / 2 \u\ A dx < (^j p\u\ 2 dx\ \\u\\ L 6 < C\\Vu\\ L 2 -)• as t ->■ oo. 

This, together with (|5.29j) and (|5.32p . proves (|1.18|) . Therefore, the proof of Theorem 11.11 is 
complete. □ 

Proof of Theorem Define the approximate data (pq V , Uq V , f &,v ) as follows: 

&,r, 35 * PO + f] S,r, . f S,n ■ , 

P = T~\ ' U =3S*Uo, f''=JS*f- 

Then, we can apply Proposition 15 . II to obtain a global smooth solution (p s,v , u 6 ' 71 ) of the Cauchy 
problem (jl.ip . (j 1 . 2 1) with the mollified data p^ ,Uq V and f S,v , which satisfies the uniform bounds 
in Proposition 13.11 Now, the remaining arguments to obtain the global weak solution and its 
asymptotic behavior are almost the same as that of [12] . The proof of Theorem 11.21 is thus 
finished. □ 
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